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Abstract: We present the results of a next-to-leading order calculation of QCD corrections 
t^J- ■ to the production of an on-shell top-anti-top quark pair in association with two flavored 

t^. ■ b-jets. Besides studying the total cross section and its scale dependence, we give several 

differential distributions. Where comparable, our results agree with a previous analysis. 
While the process under scrutiny is of major relevance for Higgs boson searches at the LHC, 
we use it to demonstrate the ability of our system built around Helac-Phegas to tackle 
complete calculations at the frontier of current studies for the LHC. On the technical side, 
we show how the virtual corrections are efficiently computed with Helac-ILoop, based on 
the OPP method and the reduction code CutTools, using reweighting and Monte Carlo 
over color configurations and polarizations. As far as the real corrections are concerned, 
we use the recently published Helac-Dipoles package. In connection with improvements 
of the latter, we give the last missing integrated dipole formulae necessary for a complete 
implementation of phase space restriction dependence in the massive dipole subtraction 
formalism. 
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1. Introduction 

Rare are processes, which received more attention than top quark pair production in various 
configurations at hadron colliders. Indeed, there is astonishing recent progress in next-to- 
leading order (NLO) [1,2] and next-to-next-to leading order (NNLO) [3-10] calculations, as 
well as next-to- next-to-leading- log resummations (NNLL) [11-13] for inclusive production. 
For this publication, however, we will mostly be interested in more exclusive channels. The 
list for the latter is just as impressive: NLO QCD corrections for the ttH signal [14-19] 
and the backgrounds from tij [20,21], tiZ [22], tibb [23,24], and most recently ttj [25]. 

The present study pertains to the tibb final state. As phenomenological motivation, 
we stress its relevance as irreducible background for light Higgs boson searches in the ttH 
channel, when the Higgs boson decays into a pair of b quarks. Indeed, realistic experimental 
analyses of this channel [26] including showering effects, b-tagging efficiencies and suitable 
cuts to reduce multijet backgrounds (among others tijj with light jets misidentified as 
b-jets), show a substantial smearing of what would be a sharp Higgs resonance peak in the 
distribution of the invariant mass of the bb system. In consequence, the knowledge of the 
backgrounds becomes crucial to claim discovery or exclusion. The dominant background is 
obviously direct production of the final state without resonances, i.e. the QCD generated 
process pp — > tibb + X. As mentioned above, the NLO computation has been completed 
only very recently in [24] . Our first goal is to confirm the results of that publication. 

The second aim is to demonstrate the power of our system based on Helac-Phegas 
[27-29], Helac-ILoop [30], CutTools [31,32] and Helac-Dipoles [33] in a realistic 
computation with six external legs and massive partons. This is the first calculation at 
this level of complexity using the OPP reduction technique [31,34] and modern unitarity 
based methods [35-41]. Indeed, computations with six external legs have only been at- 
tempted with massless partons [42-44], whereas massive partons have only appeared in a 
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more modest setting in [2]. It is fair to say, that traditional, reduction based methods, 
have proven more efficient in producing complete results until this point (see some recent 
examples [24,45-52]). 

In principle, the publication [30] has shown that the difficult virtual corrections can be 
computed within the Helac-ILoop/CutTools framework for a multitude of processes 
of practical interest, and tibb in particular. Nevertheless, it is one thing to compute a value 
for a single phase space point, and another to provide arbitrary differential distributions. 
Here we wish to convince the skeptics that the latter exercise can be mastered as well. As 
our title suggests, we are confident that we can now tackle any calculation from the "NLO 
Wishlist" [53]. 

Besides the virtual corrections, the calculation of the real radiation contributions 
is not an easy task either. Despite the announcement of several automates for Catani- 
Seymour [54,55] dipole subtraction [33,56-59], the only complete 1 and publicly available 
tool is Helac-Dipoles [33] (besides possibly Sherpa [56], which should become publicly 
available soon 2 ). In this work, we slightly increase its flexibility by allowing for phase space 
restriction [60,61] in the integrated dipoles for the general massive case (most formulae were 
known from [62,63]). 

This publication is structured as follows. In the next section, we discuss the details of 
the techniques used for the evaluation of the virtual and real corrections. Subsequently, we 
give the results for the total cross sections and differential distributions, and conclude the 
main part of the text. The Appendix contains the formulae of the massive dipole formalism 
with phase space restriction. 



2. Technicalities 

2.1 Virtual corrections 

For the calculation of virtual corrections we use Helac-ILoop [30], namely the merging 
of Helac [27-29] and the OPP [31,34] reduction code CutTools [32]. In order to further 
improve the performance and speed of the system, we make use of color, helicity and event 
sampling methods. 

The treatment of the color degrees of freedom in Helac is based on the color-connection 
representation of the amplitude. A generic QCD amplitude composed by n g gluons and n q 
quarks (and of course n q antiquarks) plus possible colorless particles can be written as 

j^n---,tn q ,ai,...,a ng 

where i and j are color indices belonging to the fundamental representation of the gauge 
group, whereas a belongs to the adjoint. Multiplying with f° ■ and summing over a, for 
each gluon, we end up with a uniform representation, namely 

Ml"-'i nq l n3 . (2.2) 



x By complete we understand phase space integration of subtracted real radiation and integrated dipoles 
in both massless and massive cases. 

2 Private communication with T. Gleisberg. 
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As it is well known, the amplitude can now be (color-)decomposed as follows 

• ••<V.M-u , (2.3) 

with k = n q + n g , and a denoting a permutation of the set {1, . . . , k}. The Feynman 
rules that allow the calculation of all A a color-stripped amplitudes in the color-connection 
representation have been described in Refs. [27-29]. The objects we are interested in, are 
the squared matrix element 

for tree order calculations, and 

for the virtual corrections, where C refers to the one-loop amplitude. 
The color sum can also be written as 

The color matrix used above is given by 

Co-,o7 = < ^<n.,ii^v a »j2 • • • < ^.,i*A CT j Ji^O ,i2 • • • <^v ,jfc ! (2-7) 
{*}.{?'} 

and is equal to 

<V = Af c m(,T '' J ' ) , (2-8) 

where m(cr, <r') count the number of common cycles of the two permutations, and Nc is 
the number of colors of the fundamental representation of the gauge group (three in our 
case) . 



Full color summation is performed in Helac by using the right-hand side of Eq. ^6 
The code generates all possible permutations of the color indices, each of them being one 
color connection. For each of the latter, using the color-connection Feynman rules, Helac 



calculates the corresponding A a color-stripped amplitude and at the end, using Eq. |2.6| , 
provides the fully color summed contribution. 

Although at tree order, practical applications are fast enough with full color summa- 
tion, at the one-loop level, one would opt for Monte-Carlo sampling over colors. This is 
done using the left-hand side of Eq. |2.6| . The idea is rather simple [64,65]: we gener- 
ate a color configuration by assigning explicit colors to the external particles. Using the 
labels 1,2,3 for the three colors in the case of QCD, a possible color assignment for a pro- 
cess like g(h,ji)g(i2,j2) -> t(i 3 ,Q)i(0,j 3 )b(u, 0)6(0, j 4 ) can be (h,...,ii) = (1,2,3,3) and 
{hi ■ ■ ■ = (3, 2, 1, 3), that is A / t 3 , 2 13 in Eq. |2.2| . Of course, color conservation requires 
that the number of times a color appears in the list of color indices (£), should be equal to 
the number of times it appears also in the list of anti-color indices (j). 
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If one now uses Eq. 2^, it is easy to see that only a few of the A a amplitudes contribute 
to a given color assignment. This means that including the color degree of freedom in our 
Monte-Carlo integration, we may reduce drastically the average number of color connections 
that are actually needed. For instance in the calculation of gg — > ttbb, the average number of 
color connections used in a Monte-Carlo sampling over color assignments, is approximately 
3 per event, resulting to almost an order of magnitude reduction in computation time with 
respect to the non sampling treatment, where all color connections, 24 in that case, have 
to be calculated per event. This is particularly important for one-loop amplitudes. It 
should be emphasized that Monte-Carlo over colors is not an approximation, it is an exact 
treatment of the color degrees of freedom. We have extensively tested that it produces the 
same results as the usual full color sum of Eq. |2.6| , see also a previous work [66] . 

The sampling over helicity is well described in Refs. [67,68]. It results to a drastic 
improvement in speed, since only one (random) helicity configuration has to be calculated 
per event, whereas in the case for instance of gg — > ttbb a full summation over helicity 
configurations will slow down the calculation by a factor that is approximately equal to 
the number of helicity configurations, namely 32 in that case. 

Finally the actual calculation of the virtual corrections is organized using a re- weighting 
technique [69,70]. To explain how this works, let us start with the following equation 

°a° +V = f dxidx 2 d^ m f a ( Xl )f b (x 2 ) {\M\ 2 + MC + M*C) , (2.9) 

which gives the sum of leading order (LO) and virtual (V) contributions for a scattering 
ab — > m-particles. It can be re-written as 

// SAC* A- SA* C\ 
dx l dx 2 d^ m f a (x 1 )f b {x 2 )\M\ 2 (l + -±p J . (2.10) 

Since C is a time consuming function one would like to calculate it as few times as possible. 
To this end a sample of un-weighted events is produced based on the tree order distribution, 
namely 

1 da LO 

g(X) = g( Xl ,x 2 , d> ro ) = ab , (2.11) 

a^ u dx\dx 2 d® m 

satisfying J dXg(X) = 1. The sample S of un-weighted events has the following property, 

0(Xi) = [ dXg{X)0{X) , (2.12) 
5 ies J 

where the equality should be understood in the statistical sense, and O(X) is any well- 
defined function over the integration space. Now it is trivial to see that if 



then 



1 ct lo+v 



Ns 7ts °oh 
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In practice the sample of tree order un-weighted events includes all information on the 
integration space, namely, the color assignment, the (random) helicity configuration, the 
fractions x\ and X2 and the m— body phase-space. For future convenience it is produced 
in a standard Les Houches format [71]. One-loop contributions are only calculated for this 
sample of un-weighted events, and the weight assigned to each of those events is given by 

MC* + M*C 

The total virtual contribution can now be easily estimated by 

a v = (w) a B , (2.16) 

where a B is the born cross section, already included in the Les Houches file. Moreover, the 
sample of events including the information on w, can be used to produce any kinematical 



« = VT7T2 • ( 2 - 15 ) 



distribution, according to the Eq. 2.12. 



In our application, the speed-up factor obtained in this way varies between 10 2 and 10 3 . 
In fact, an event in the gluon-gluon channel costs about a second on a 3 GHz machine, and 
few permille level cross sections are obtained by re-weighting about fifty thousand events. 
The smooth distributions of the next sections are obtained on samples of two hundred 
thousand events. We stress it for the non-expert reader that all timings are subject to 
drastic future improvement, but it is easy to see that the complete calculation of the 
virtual corrections is a matter of one or two days on a single machine. We consider this a 
practical proof of the power of our approach. 

Finally, for the expert reader this time, we stress that the numerical stability of the 
virtual correction evaluation is checked by performing a gauge independence test on each 
event. 

2.2 Real radiation 

The real radiation contribution to the process is obtained using the Catani-Seyour dipole 
subtraction method [54] in the massive version as described in [55] and extended for ar- 
bitrary polarizations in [33]. We use a phase space restriction on the contribution of the 
dipoles as originally proposed in [60,61]. Most of the formulae needed beyond that work 
to account for massive partons have been presented in [62,63]. In fact, the only missing 
integrated dipoles correspond to the final emitter and final spectator case, when both are 
massive. This situation occurs in our calculation, when a top quark emits a soft gluon, 
which is then absorbed by the anti-top-quark acting as the spectator of the Catani-Seymour 
formalism. As advertised in the introduction, we give the complete set of expressions in 
the most general case in the Appendix. Let us stress at this point, that, similarly to most 
authors, we do not use finite dipoles regularizing the quasi-collinear divergence induced by 
both top quarks moving in the same direction, even though they are implemented in the 
software. Due to the large top quark mass, they are not needed to improve the numerical 
convergence. 

At this point let us remind the reader that the phase space restriction on the dipole 
phase space is defined differently depending on whether the spectator and emitter are in 
the final or initial states. To be more specific a given dipole contributes as follows 
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1. for final- final dipoles, if 



Vij,k < C , (2-17) 

where 

ViJ,k = > (2- 18 ) 

PiPj + PiPk + PjPk 

and Pi, pj and pk are the momenta of the emitter pair and of the spectator respectively. 

2. for final- initial dipoles, if 

1 — x ij,a < a max > (2-19) 

where 

_ PaPi + PaPj ~ PiPj + 5 ("4 ~ m i ~ m j) (22Q) 
PaPi + PaPj 

and Pi,Pj and p a are the momenta of the emitter pair and of the spectator respectively. 
rriij is the emitter particle mass. 

3. for initial-final dipoles, if 

m < aZ x , (2.21) 

where 

Ut = ggg , (2.22) 

PiPa + PjPa 

and p a , Pi are the emitter pair momenta, with p a in the initial state, and pj is the 
final state spectator momentum. 

4. for initial-initial dipoles, if 

^ - rf < - (2-23) 

where p a ,Pi are the emitter pair momenta, with p a in the initial state, and pb is the 
initial state spectator momentum. 

In our implementation of the subtracted real radiation the a parameters can be varied 
independently of each other. On the other hand, we slightly simplified the implementation 
of the integrated dipoles by assuming a^cL = a max = a lnax = a max = a max- We consider 
two extreme choices, namely a max = 1 and a max = 0.01. The first one corresponds of 
course to the original formulation of [55]. 

Even though the use of a phase space restriction in the dipole formalism is rather 
widespread currently, we wish to make a few remarks on its advantages and disadvantages 
in the present setting. On the side of advantages there are in fact three reasons 
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Figure 1: Distribution of the invariant mass m b i of the bottom- anti-bottom pair for pp — > ttbb+X 
at the LHC for different parts of the real radiation contribution with different choices of a max , 
Oimax = 1 o,nd a max = 0.01. The red solid line corresponds to the sum of all contributions, the 
blue dashed line represents the dipole subtracted real emission, the cyan dot- dashed line corresponds 
to the sum of the K and P insertion operators, and finally the green dotted line represents the I 
insertion operator. The sum of all the contributions for the two different choices of a max is depicted 
below. 



1. Our phase space generator, Phegas [72], uses multi-channel optimization [73], where 
the phase space density of a given channel corresponds to a product of Feynman 
diagram denominators elevated to some arbitrary power, which turned out to be 
best chosen relatively close to unity. This clearly reproduces the peaked behavior in 
the collinear and soft limits, with the integrable square root singularity after dipole 
subtraction. On the other hand, the distortion of the remaining (non-singular, but 
nevertheless present) peaks in the dipole subtraction terms due to phase space remap- 
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ping is not taken into account at all. Therefore, restricting the kinematics to regions 
very close to the singular limits keeps the distortion small and allows to obtain a 
maximum gain from the available channels. 

2. In close relation to the previous issue, we observe that there is more than a factor of 
three less events accepted if the phase space restriction parameter a max is equal to 
0.01. The difference is due to the fact, that after phase space remapping, an event 
which would not be accepted by real radiation cuts, may pass the cuts of the dipole 
jet function if a max = 1. This phenomenon is called missed-binning. It is important 
to note that despite more accepted events, the absolute error of the final result is 
only slightly better than for a max = 0.01. In fact, to obtain the same absolute error, 
about twice less events are needed with the latter choice. This is is one of the two 
speed-up factors. 

3. On the average much less dipole subtraction terms are needed per event with a max = 
0.01, since the collinear limit singles out a pair of partons, and the soft limit requires 
the sum over all dipoles involving the soft parton as emitter only. This constitutes 
the second speed-up factor. 

Unfortunately, having a max = 0.01 implies large cancellations between the dipole 
subtracted real radiation and the integrated dipole contribution. With this choice of a max , 
the former overshoots the complete result by a factor of almost three. It is probably safe 
to state that a slightly larger value of a max , such that the real radiation contribution 
(subtracted with dipoles) would be similar in value to the final result for the sum of real 
radiation and integrated dipoles, would be more advantageous. We leave this issue to 
future studies. 

It is interesting to observe the different contributions on a chosen distribution. We 
illustrate such a decomposition in Fig. [l], for the invariant mass of the bb system. All 
parameters and cuts are defined in the next section. We separate the contribution of the I 
operator and of the sum of the K and P operators defined in [55]. When introducing the 
otmax parameter dependence, we have again some freedom, and we choose our operators 
such, that the dependence of the I operator is exactly as in the massless case of [61], except 
for the final-final integrated dipoles for which we give the formulae in the Appendix. Clearly, 
the mentioned large cancellations for a max = 0.01 involve all three contributions, but the 
integrated dipoles prove to have much better convergence and the final error is entirely 
dominated by the statistically demanding real radiation. The plot corresponding to the 
sum of all contributions, also presented in Fig. |], proves the agreement between the results 
for the two parameter choices. We will show the agreement on the total cross sections in 
the next section. 

3. Results 

We consider the process pp — > ttbb + X at the LHC, i.e. for y/s = 14 TeV. For the top- 
quark mass we take mt = 172.6 GeV, whereas all other QCD partons including b quarks 
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are treated as massless. Mass renormalization is performed in the on-shell scheme. All 
final-state b quarks and gluons with pseudorapidity |?7| < 5 are recombined into jets with 
separation y Acf) 2 + Ay 2 > D = 0.8 in the rapidity azimuthal-angle plane via the IR-safe 
/^-algorithm [74-76] . Moreover, we impose the following additional cuts on the transverse 
momenta and the rapidity of two recombined b-jets: pp,b > 20 GeV, \yb\ < 2.5. The 
outgoing (anti)top quarks are neither affected by the jet algorithm nor by phase-space 
cuts. The separation between the b-jets, AR bb > D, implied by the jet algorithm, together 
with the requirement of having both b-jets with px b > 20 GeV sets an effective lower limit 
on the invariant mass m bb 

(m bb ) min = (pT,b)minV 2 0- ~ cosD) 15.6 GeV , (3.1) 

which screens off the collinear singularity. This is the reason, why we don't need any dipoles 
for the gluon splitting into a bb pair. 

We consistently use the CTEQ6 set of parton distribution functions (PDFs) [77, 78] , 
i.e. we take CTEQ6L1 PDFs with a 1-loop running a s in LO and CTEQ6M PDFs with 
a 2-loop running a s in NLO, but the suppressed contribution from b quarks in the initial 
state has been neglected. The number of active flavors is Np = 5, and the respective 
QCD parameters are Ag = 165 MeV and A5 = 226 MeV. In the renormalization of the 
strong coupling constant, the top-quark loop in the gluon self-energy is subtracted at zero 
momentum. In this scheme the running of a s is generated solely by the contributions of 
the light-quark and gluon loops. By default, we set the renormalization and factorization 
scales, fj,R and fj,p, to the common value (jlq = rrif- 

We would like to stress that the above parameters correspond exactly to those assumed 
in the analysis of [23, 24], and are essentially based on [26]. It is clear that there are 
many interesting phenomenological analyses that can be performed using our system with 
different cuts, but, as explained in the introduction, our main goal is to demonstrate its 
correctness and efficiency. To this end we want to be able to compare with the previous 
study. 

We begin our presentation of the final results of our analysis with a discussion of the 
total cross section at the central value of the scale, fiR = \ip = mt- The respective numbers 
are presented in Tab. [l] for the two choices of the a max parameter. We also single out the 
quark channel (although its contribution beyond leading order is negligible for any practical 
study) because we can compare our results with [23] . Clearly, we observe perfect agreement 
within statistical errors between all independent evaluations. Notice that we quote smaller 
statistical errors than [24] for the complete proton-proton scattering cross section, mostly 
because we performed many experimental computations of the different sub-parts and the 
accumulated statistical sample is fairly sizable. At the central scale value, the full cross 
section receives a very large NLO correction of the order of 77% which is mainly due to 
the gluonic initial state as stressed previously in [24]. 

Subsequently, we turn our attention to the scale dependence, which is given in Tab. || for 
the total cross section for pp — ► tibb + X at the LHC at LO and NLO with //r = fip = £mt 
for a few distinct values of £. As expected, we observe a reduction of the scale uncertainty 
while going from LO to NLO. Varying the scale up and down by a factor 2 changes the 



-9- 



Process 


4°3,24] ^ 


CT L0 [fb] 


rr NLO ffbl 
CT [23,24] ^ 


aNLO_ [fb] 

"max — J- L J 


a Qmi°=o.oi [fb] 


qq — ► 


85.522(26) 


85.489(46) 


87.698(56) 


87.545(91) 


87.581(134) 


pp — > 


1488.8(1.2) 


1489.2(0.9) 


2638(6) 


2642(3) 


2636(3) 



Table 1: Cross sections for pp — > ttfefe + 1 d i/ie _Li/C at LO and NLO for the scale choice 
Hf = Hr — rat, in comparison with the results of Refs. [23,24]- The statistical errors are quoted in 
parentheses. 





1/8 • mt 


1/2 • rn t 


1 • m t 


2 • m t 


8 • m t 


a LO [fb] 


8885(36) 


2526(10) 


1489.2(0.9) 


923.4(3.8) 


388.8(1.4) 


a NLO [fb] 


4213(65) 


3498(11) 


2636(3) 


1933.0(3.8) 


1044.7(1.7) 



Table 2: Scale dependence of the total cross section for pp — > ttbb + X at the LHC at LO and 
NLO with \iR — [if = £mt for a few distinct values of £. 



cross section by +70% and -38% in the LO case, while in the NLO case we have obtained 
a variation of the order +33% and -27%. Our findings can be summarized as follows 

o£g(LHC, mt = 176.2 GeV, CTEQ6L1) = 1489.2 ™ fb , (3.2) 

a^ L 5 °(LHC, m t = 176.2 GeV, CTEQ6M) = 2636 ^ Sj fb . (3.3) 

As can been easily seen from Tab. |2], choosing a scale too small (too large) will corre- 
spond to a very large negative (positive) correction. Unfortunately, this large scale variation 
and the size of the corrections themselves, imply that if a meaningful analysis were required 
in the present setup, a full NNLO study would be indispensable. As the latter will remain 
out of reach in the nearest future, it seems that, as already suggested in [24], additional 
cuts must be introduced in order to reduce the NLO corrections. Only then will this pro- 
cess, which constitutes the main irreducible background, not put in danger the feasibility 
of Higgs boson searches in the ttH channel. 

In Fig. [2] we show the result for the scale dependence graphically. Since we confirm the 
findings of [23,24] for all other numbers, we will not present the scale dependence with the 
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Figure 2: Scale dependence of the total cross section for pp — ► tibb + X at the LHC with fiR 
III £to 4 . On the upper panel, the blue dashed curve corresponds to the leading order, whereas the 
red solid one to the next-to-leading order result. The lower panel shows the scale dependence of the 
next-to-leading order cross section (red solid curve) decomposed into the contribution of the virtual 
corrections (green dashed curved) and real radiation (cyan dash-dotted curve). 



renormalization and factorization scales varied independently, which can be found in that 
work. On the other hand it is entertaining to see the final scale dependence for both scales 
equal emerge out of the two contributions (virtual and real), as also depicted in Fig. |2[ 
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Figure 3: Distribution of the invariant mass m b g of the bottom- anti-bottom pair (a), distribution 
in the transverse momentum px bl of the bottom- anti-bottom pair (b), distribution in the rapidity y b i 
of the bottom- anti-bottom pair (c) and distribution in the transverse momentum px b of the bottom 
quark (d) for pp — » ttbb + X at the LHC at LO (blue dashed line) and NLO (red solid line). All 
distributions have been obtained with a max = 0.01. 

Of course, the separation is entirely unphysical, but well defined once we state that we 
use what is now called the 't Hooft-Veltman [79] version of the dimensional regularization, 
with the integrals as defined in [30]. The large cancellation for small scale values is the 
source of the rising integration errors quoted in Tab. |2[ 

The entirely new result, which has not been presented in the literature until now, 
are the differential distributions for the four simplest observables, namely the invariant 
mass, transverse momentum and rapidity of the two-6-jet system, as well as the transverse 
momentum of the single 6-jet. These results can be found in Fig. || and have been obtained 
with a max = 0.01. The histograms for a max = 1 are of similar quality, but we have 
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m bb [GeV] y bb 

Figure 4: The ratio of the NLO and LO distributions for pp — ► tibb + X at the LHC for the 
invariant mass m hb (left panel) and the rapidity y b i (right panel) of the bottom- anti-bottom pair. 



refrained from averaging over the two statistically independent evaluations. Clearly, the 
distributions show the same large corrections, which turn out to be relatively constant 
contrary to the case of quark initial states as shown in [23]. 

The histograms can also be turned into dynamical K-factors, which we show in two 
cases in Fig. f|, namely for the invariant mass m bb and the rapidity y bb of the two-6-jet 
system. For those two cases, they are simply defined as 

_ da^°/dm bb 

K ^b)- da LO/ dmbl > W 



and 



R( , _ da^Q/dy bb 

K ^b)- d(J LO/ dybl > W 



respectively. As already anticipated above, we notice that they have a relatively small 
variation, when compared with their size. 



4. Conclusions 

We have presented a complete study of the pp — > tibb + X process at the LHC at the 
next-to leading order in QCD. Our results agree with the previous study of [23,24], for all 
the numbers given in that work. This fact alone constitutes a powerful proof that our fully 
automatic system can be applied in analyses of realistic processes for the Large Hadron 
Collider (and for any other collider for that matter). 

One conclusion, which follows from the large observed corrections, is that in view of 
the importance of the process under study as background for Higgs boson production in 
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association with a top quark pair, a more detailed phenomenological analysis is necessary. 
We postpone such work for the future. 

As a completely technical detail, our work provides the last missing formulae for the 
dipole subtraction formalism with a dependence on a dipole phase space restriction param- 
eter. This allowed for internal tests, as well as useful speed-ups of the calculation. 



Acknowledgments 

The work of M.C. was supported by the Heisenberg Programme of the Deutsche Forschungs- 
gemeinschaft. G.B., C.G.P., R.P. and M.W. were funded in part by the RTN European 
Programme MRTN-CT-2006-035505 HEPTOOLS - Tools and Precision Calculations for 
Physics Discoveries at Colliders. M.W. was additionally supported by the Initiative and 
Networking Fund of the Helmholtz Association, contract HA-101 ("Physics at the Teras- 
cale"). R.P. thanks the financial support of the MEC project FPA2008-02984. 



A. Integrated dipoles with emitter and spectator in the final state 

In this Appendix, we give the formulae necessary for a complete implementation of the 
phase space restriction in the dipole subtraction formalism. We follow closely the notation 
of [55], and present only the minimum necessary for an implementation in a numerical 
program. 

As stated before, the only missing ingredient concerns the case with a final state emitter 
and final state spectator. The relevant kinematical variable is 

Vij,k = ; — , A.l 

PiPj + PiPk + PjPk 

where Pi,Pj are the momenta of the emitter pair and p^ is the momentum of the spectator. 
We consider the case, where mf = pj = 0, and thus ignore the splitting of a gluon into a 
pair of massive quarks. Note that the splitting of a heavy quark into a heavy quark and a 
gluon in the presence of a massless spectator, as well as the splitting of a massless quark 
in the presence of a massive spectator have already been covered in [63]. The remaining 
options can be found below. 

The upper limit on y^^ without a phase space restriction is 

1 - c j - n 

where fij t k = mj,k/ \f& , with Q = Pi + pj + Pk- The restriction is imposed by adding the 
following condition 

Vij,k < a . (A.3) 
Clearly, we will only observe a modification if 

a<y+ . (A.4) 



-14- 



In such a case, we write the integrated dipole function corresponding to Eq. 5.22 of [55] as 

hj,k{ e ) A?',fc( e > a ) = Iij,k( e ) + ^hj,k{a) ■ (A.5) 

Notice that, since setting a minimum value of y^^ screens from all divergences, the addi- 
tional term A 7^ (a) is finite and can be obtained by integration over the four-dimensional 
phase space, thus simplifying substantially the calculation. Indeed, we will use the phase 
space of Eq. 5.11 from [55], with e = and the following insertion 

J[dpi(Pij,Pk)] -> J[dPi(Pij,Pk)]0(yij,k> a) . (A.6) 
We first give the result for the eikonal integral defined as 

with Zj = PjPk/(PiPk +PjPk)- It can most conveniently be obtained by a variable change, 
and we therefore parametrize the formulae with 



x = y + - a + 



1 4:fj?- fjP 



We did not attempt to minimize the number of dilogarithms and our result contains sixteen 
of them 

A/- k = J-f-Li 2 ( + Li2 ( °) + Li2 ( _ Li2 ( *+ 



Vjk\ \a + x + J \a + x+J \x + — bj \x + — b 

f c + x\ f c \ T . /x_ — x\ ( X- 

+Li 2 - Li 2 + Li 2 - Ll 2 



c + x + J \c + x + J \x_+a/ \x_ + a 

,. , b—x\ r . ( b \ T . ( X- — x\ T . / x. 
-Li 2 + Li 2 - Li 2 + Li 2 





— X 


X+ 


-b 


x_ 


— X 


x_ 


+ a 


x_ 


— X 



b — x-J \b — x-J \x- + cj \x- + c 

+Li 2 ( 7 — - ) - Li 2 It—) ~ Li 2 ( < —^- I + Li 2 



b + aj \b + a J \c — aj \c — a 

, \ i /(a — c)(x+ — x)\ , /x, / ( a - c ) a:; + 

+ log (c + x) log 7 T - lo §( c ) lo S 



a + x)(c + x+) / \a(c + x + ) 

.'(a + x)(x_ (a{x--b) 

+ log (6 - Xj log — r - log(O) log 1 



(a + fr)(x_ — x) J \{a + b)x 

- log ((a + x)(6 - x+)) log (x+ — x) + log (a(6 - log(x+) 

, /a, / (a + x)x + x_ \ , /x_ — x\, fc + x 
+ log(d) log — i (-t + log log 



where 



a(x+ — x)(x_ — x) / \ x_ / ya + x-^ 

+ 2~ l0g (~~a~~) bg + x ^ a + x +) 2 )) ' ( A - 9 ) 
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b = 



2(1 - fi k ) 
1-^ 



2 ' 

4 



c = 



2/i fc (l - Hk) 
^(1-M?-Mfc) . 



(A.ll) 
(A.12) 
(A.13) 



and 



(1-^0 -^±a/a(i,/4^) 



(A.14) 
(A.15) 



with A the Kallen function 



A(x, y, z) = x 2 + y 2 + z 2 — 2xy — 2xz — 2yz 



(A.16) 



Besides the eikonal integral we will also need the collinear integrals (for this and fu- 
ture works with massive quarks), which are implicitly defined as (in complete analogy to 
Eqs. 5.23-5.25 of [55]) 



2A/ cik + AI$} k 



AI, 



qq,k 



TrAP 



coll 



AI, 



2C A 



qq,k J 

2A/ eik + A/ 



coll 

gg,k 



(A.17) 
(A.18) 
(A.19) 



The respective results read refers to the mass of the heavy quark) 



AI, 



gQ,k 



AI, 



qq,k 



AI, 



gg,k 



2 (1 + a) + T^ 



2(2 - 2fi 2 - fi k ) 



+ 



(1 — a)// 



1 - M 2 _ ^ 2(jJ 2 + a(1 _ ^2 _ ^2)) 



(A.20) 



" 21 °g - -^5 —\+— -9. — log 1 



(1 - Mfc ) 2 - ^ 2 I 2(1 - ^ 2 - fi) 



2 ( I- fjL k - a(l + /i fc ) 1q / a(l + ji k ) 

3 \ 1 + ^fe V 1 -^ 



(1-W) 2 



(A.21) 



+ K 



2\ 2^ 2 

- log 



(l-a)(l + /i fc ) 



37 i-^i 

11 (\- n k - a(l + /U fc ) 



2/Ufc 



6 



+ log 



«(1 + Atfc) 
1 - Mfc 



3/ 1 



1 + Mfc 

^ w A 1 -«)(!+ A*fc) 

2 1U 6 



(A.22) 



2/Xfe 
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